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=  Area  of  the  Adsorbent 
=  Langmuir  isotherm  fitting  parameter  (1/MPa) 

=  Cumulative  (integral)  energy  distribution 
=  Single  probability  density  function 
=  Joint  probability  density  function 
=  Exponential  curve  fitting  parameters 
=  Freundlich  parameters 
=  a^la^ 

=  Global  amount  adsorbed  (moles/1^) 

=  Amount  adsorbed  in  local  homogoieous  patch  (moles/1^) 
=  Partial  pressure  (MPa) 

=  Saturation  adsorption  value  (moles/1^) 

=  Universal  gas  constant  (kJ/mole»K  or  MPa»L/mole*K) 
=  Constant  separation  factor  isotherm  parameter 
=  Temperature  (K) 

=  Adsorbed  phase  mole  fraction 
=  Vapor  phase  mole  fraction 
=  Total  number  of  components  in  the  mixture 


GREEK  LETTERS 

Yj  =  Activity  coefficient  in  the  adsorbed  phase 
€  =  Adsorption  site  energy  (U/mole) 

/tp  =  p*  cratral  moment  of  a  probability  doisity  function 
n  =  Spreading  pressure 

p  =  Site  matching  parameter  or  covariance  of  the  probability  density  function 

o  =  Square  root  of  the  variance  of  a  probability  density  function 

o*  =  Excess  area  of  mixing  (kg/mole) 

4>  =  Fugacity  coefficient 

ijr  =  Spreading  pressure  function  nA^RT  (moles/kg) 
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LIMITATIONS  IN  THE  APPLICATION  OF  THE  ADSORBED 

SOLUTION  THEORY 


INTRODUCTION 

The  need  exists  for  a  method  to  accuratdy  model  mixture  adsorption  equilibria.  This 
need  arises  not  only  from  an  academic  desire  to  understand  the  fundamentals  of  the  adsorption 
process,  but  also  from  the  more  pragmatic  requirement  to  provide  the  nucleus  of  mathematical 
models  which  describe  dynamic  adsorption  systems.  In  the  course  of  generating  results  for  a 
dynamic  adsorption  system  using  accqjted  mathematical  models,  adsorption  isotherm 
relationships  are  generated  hundreds  of  thousands  of  times.  This  points  to  the  importance  of 
gmoating  equilibrium  adsorption  data  with  a  minimum  of  computational  burdra  and  a  high 
degree  of  accuracy. 

Several  models  have  bewr  developed  to  describe  adsorption  systems.  The  simpler 
models  often  are  limited  in  their  applicability  and  accuracy,  while  the  more  complicated 
systems  become  too  cumbersome  or  unsolvable  under  some  conditions.  Several  of  the  models 
are  ultimately  correlative  in  nature.  While  these  models  provide  important  information  in 
understanding  the  mixture  adsorption  process,  their  utility  is  limited  due  to  the  paucity  of,  and 
difficulty  in  obtaining,  mixture  adsorption  data.  As  a  result,  models  that  can  predict  mixture 
adsorption  from  single  component  adsorption  data  are  becoming  more  popular. 

This  study  will  investigate  one  of  the  mote  popular  mixture  adsorption  theories,  the 
Adsorbed  Solution  Theory  (AST)  (Myers  and  Prausnitz,  1965),  and  its  ability  to  predict 
mixture  adsorption  from  single  component  data.  In  an  effort  to  isolate  pure  theory  from  the 
effects  and  inherent  errors  associated  with  experimental  data,  this  project  was  divided  into  two 
tasks:  A)  The  development  of  a  computationally  fast  adsorption  isotherm  model  that  can  be 
2q>plied  to  systems  whidi  demonstrate  energetic  heterogeneity;  B)  Limitations  in  the  application 
of  the  adsorbed  solution  theory.  NRL  Memorandum  Report  6116-97.7946  (Nfrituszko  and 
Lamontagne,  1997)  detailed  the  development  of  part  A.  This  report  (part  B)  will  cover  the 
complications  and  problems  encountered  wh^  applying  the  AST  to  systems  that  have 
unfavorably  adsorbed  components.  Since  this  work  focuses  primarily  on  the  AST,  a  short 
background  section  of  this  theory  is  provided.  For  an  in-depth  explanation  of  AST  as  applied 
to  this  problem,  the  reader  is  referred  to  Matuszko  and  Lamontagne  (1997). 

ADSORBED  SOLUTION  THEORY 

BACKGROUND 

The  Adsorbed  Solution  Theory  (AST)  is  the  most  widely  used  model  to  describe 
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adsorption  of  gas  phase  mixtures.  The  AST  is  a  translation  of  three  dimensional  solution 
thermodynamics  to  the  two  dimensions  of  an  adsorbed  surface.  Three  dimensional  properties 
such  as  pressure  and  volume  become  the  two  dimensional  properties,  spreading  pressure  and 
area.  In  addition,  mixing  properties  associated  with  nonidealities,  such  as  activity  coefficients 
and  excess  volume,  become  adsorbed  phase  activity  coefficients  and  excess  area.  The  popular 
theories  for  incorporating  nonidealties  into  the  adsorption  system  detailed  in  Matuszko  and 
Lamontagne  (1997),  assume  that,  while  the  overall  mixture  interactions  may  be  affected  by  the 
adsorbent  surface,  the  relative  interactions  between  adsorbates  remain  the  same  as  those  found 
in  a  liquid  mixture.  It  is  expected,  therefore,  that  adsorbed  mixtures  should  behave  similarly 
to  the  analogous  liquid  mixtures. 

The  adsorption  of  water  and  organic  mixtures  is  a  case  where  the  translation  from 
solution  properties  to  adsorption  properties  is  of  particular  interest.  As  liquids,  water  and 
organic  mixtures  are  so  highly  non-ideal  that  two  sqiarate  phases  are  formed.  The  translation 
of  extremely  non-ideal  behavior  to  the  adsorbed  phase  has  not  provided  adequate  results 
(Matuszko,  1992),  and  therefore  requires  a  more  thorough  investigation. 

To  simplify  the  investigation  of  nonidealities,  Talu  and  Zwiebel  (1986)  developed  a 
method  which  strips  all  but  the  most  basic  assumptions  governing  adsorbate/adsorbate 
interactions  fi-om  the  AST.  The  only  inherent  assumption  in  the  method  is  that  the  adsorption 
system  can  be  described  with  spreading  pressure  dependent  activity  coefficients.  Otherwise, 
all  assumptions  arising  from  the  choice  of  a  theory  for  describing  excess  Gibbs  free  energy, 
and  hence  activity  coefficients  and  excess  area  of  mixing,  are  removed.  This  method  is  a 
powerful  tool  in  elucidating  the  exact  relationships  between  adsorbates.  The  method  involves 
the  direct  calculation  of  mixture  spreading  pressures  from  binary  data  and  uses  them  to 
generate  adsorbed  phase  activity  coefficients  and  excess  areas  of  mixing.  These  mixing 
parameters  are  then  used  to  predict  the  performance  of  ternary  or  other  systems  which  have  no 
mixture  data. 

This  analysis  uses  Talu  and  Zwiebel's  (1986)  method  to  generate  adsorbate/adsorbate 
interaction  paramet^  for  a  mixture  of  perfluorocyclobutane  (PFCB)  and  water  on  BPL 
activated  carbon.  The  purpose  of  this  effort  is  to  determine  the  effect  the  adsorption  process 
has  on  the  overall  value  of  the  activity  coefficients  and  excess  area  of  mixing.  While  it  is 
accepted  that  the  relative  values  of  these  parameters  remain  the  same  as  those  found  in  a  liquid 
mixture,  the  overall  impact  of  the  adsorption  process  remains  unknown.  In  addition,  as  a 
byproduct  of  the  results  of  this  analysis,  a  mathematical  investigation  is  made  into  the 
application  of  the  AST  to  an  organic/water  adsorption  system. 

ANALYSIS  OF  PRIOR  EXPERIMENT 

Adsorption  isotherm  data  for  perfluorocholorobutane  (PFCB)  and  water  on  BPL 
activated  carbon  at  room  temperature  were  obtained  from  Matuszko  and  LeVan  (1994).  The 
apparatus  and  procedure  used  for  the  experiments  is  described  in  Farris  et  al.  (1993).  The 
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data  were  taken  at  a  constant  vapor  phase  mole  fraction  so  that  experimental  mixture  spreading 
pressures  could  be  obtained  by  direct  integration  of  the  data.  In  addition,  the  experiments 
were  controlled  so  that  each  point  would  fall  on  the  adsorption  leg  of  the  water  adsorption 
hysteresis  loop.  These  data  can  be  found  in  Table  1. 

According  to  Talu  and  Zwiebel  (1986),  when  data  are  taken  in  the  constant  vapor  phase 
mole  fraction  plane,  the  experimental  spreading  pressure  (i|r)  can  be  obtained  from: 

^  ~  io  (1) 


This  integral  is  evaluated  by  fitting  the  mixture  data  points  from  Table  1  to  the  curve: 
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Where:  H  =  0.12866 

L  =  2087.2 


Equation  2  is  plugged  into  Equation  1,  integrated,  and  the  spreading  pressure  for  each  mixture 
data  point  calculated.  Once  the  experimental  spreading  pressure  (ijr)  is  obtained  for  each  point, 
the  AST  reference  values  (nj®  and  Pj”)  can  be  obtained  using: 
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ITie  integral  in  Equation  3  is  solved  for  n  ®  using  Equation  4  (the  single  component  isotherm 
correlation  for  each  specie).  For  the  purposes  of  this  effort,  a  three  parameter  modified 
Antoine  equation  (Hacskaylo  and  LeVan,  1985)  was  used  for  PFCB  and  a  four  parameter 
polynomial  was  used  for  water.  The  forms  for  these  equations  are  identical  to  those  used  in 
Matuszko  (1992).  The  corresponding  P®  value  is  obtained  from  Equation  4.  The  reference 
values  are  combined  with  the  experimental  mixture  values  using: 


Pi  =  X,Phi 


Y,  = 


xp; 


(5) 


3 


n 


+  a' 


tot 


'=l  fli 


a*  =  -i-  -  E 


2  T 


-  «  ^ 

.=1  n, 


(6) 


The  values  calculated  for  o®  and  Yj  can  be  found  in  Table  1. 


TABLE  1:  PFCB  and  Water  mixture  adsorption  data  and  AST  calculatedj^^ 


p. 

p. 

n. 

_ii, — 

x._ 

- !l!_ 

o* 

- - 

- Yi - 

4.613C-04 

1.016C-05 

0.300 

0.176 

0.6296 

0.1660 

^.6752 

9.6013 

11369.8 

7.792c4)4 

1.530c^5 

1,167 

0.250 

0.8236 

03480 

-2.5326 

5.5778 

13909.6 

1.004C-03 

2.053c4)5 

0.719 

0.256 

0.7372 

0.5331 

-13005 

5.1050 

6667.0 

1.300C-03 

2.498e-05 

2.000 

0.305 

0.8677 

0.8852 

4).8345 

3.3254 

6810.5 

1.449C-03 

2.990C-05 

1.651 

0.292 

0.8498 

1.1362 

4).4877 

2.9499 

4443.8 

1.680C-03 

3.557C-05 

3,915 

0.310 

0.9267 

1.6541 

-0.3836 

2.1884 

4833.6 

1.942C-03 

4.113c^5 

8.022 

0.194 

0.9763 

2.5251 

-0.2162 

1.6597 

6014.6 

2.428C-03 

4.982C-05 

22.46 

0.020 

0.9991 

5.5864 

-0.0254 

1.2897 

15322.6 

1  =  Water  2  =  PFCB 


The  most  important  information  found  in  Table  1  are  the  calculated  o®  values.  These 
values  represent  the  change  in  adsorption  area  occupied  by  the  components  as  a  result  of  their 
interaction  with  the  other  specie.  Since  in  an  ideal  mixture  there  is  no  area  change  on  mixing, 
o'  is  always  zero.  The  table  shows,  however,  that  o'  is  always  negative  for  PFCB  and  water 
on  BPL  carbon.  Negative  o'  values  indicate  that  the  mixture  is  occupying  less  adsorption 
space  than  would  an  ideal  solution.  This  result  would  indicate  that  PFCB  and  water  associate 
and  attract  each  other  in  the  adsorbed  phase  and  therefore  occupy  less  adsorption  space. 
Considering  that,  in  the  liquid  phase,  PFCB  and  water  are  highly  non-ideal  and  form  two 
separate  phases,  the  results  above  seem  incorrect  and  require  further  investigation. 

THEORY 

To  address  the  excess  area  anomaly  found  with  PFCB  and  water,  a  more  thorough 
investigation  into  the  fundamentals  of  the  AST  was  initiated.  An  approach  was  developed  to 
examine  the  implications,  in  general,  of  using  the  AST  on  a  mixture  which  displayed 
unfavorable  single  component  adsorption  behavior.  Both  Freundlich  and  constant  separation 
factor  isotherms  are  considered  below. 
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Freundiich  Isotherm 


Because  of  its  mathematical  simplicity,  and  its  ability  to  easily  represent  both  favorable 
and  unfavorable  adsorption  isotherms,  the  Freundiich  isotherm  was  selected  for  this  model 
study.  The  Freundiich  isotherm  takes  the  form: 


n,  = 


Vu, 


(7) 


When  Uj  is  greater  than  1,  the  isotherm  is  favorable,  and  when  U;  is  less  than  1,  the  isotherm  is 
unfavorable.  To  further  simplify  the  calculations,  no  mixing  nonidealities  will  be  assumed 
(i.e.  the  Ideal  Adsorbed  Solution  Theory  --  lAST  will  be  used). 

The  lAST  calculations  for  a  binary  system  have  been  performed  by  Crittenden  et  al. 
(1985)  on  a  Freundiich  isotherm  taking  the  Equation  7  form.  Crittenden  et  al. ,  after  much 
algebra,  produced  the  following  lAST  relationship: 


If  a  derivative  with  respect  to  nj  is  taken  the  equation  becomes: 
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For  normal,  ideal,  binary  isotherm  behavior,  P2  must  increase  if  n,  increases  holding  n2 
constant.  In  an  ideal  adsorbed  solution,  the  adsorbate-adsorbate  interactions  should  be  limited 
to  the  occupation  of  adsorption  space.  Any  increase  in  the  amount  adsorbed  of  one  component 
takes  up  room  and  tends  to  crowd  out  the  second  component.  For  the  adsorbed  amount  of  the 
second  component  to  remain  unchanged,  there  should  be  an  increase  in  the  component  2 
pressure.  Therefore,  the  quantity  in  the  brackets  in  Equation  10  should  always  greater  than 
zero: 


1  w.w,  u, 

"i-(— - )( - +"2)  =  “i-(— ^1+^2)  >0 
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(11) 
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OR 


(12) 


If  u,  and  U2  are  both  greater  than  1,  Equation  12  is  always  satisfied.  If,  however,  Ui  and/or  Uj 
are  less  than  1,  then  die  condition  may  not  be  satisfied.  If  the  condition  is  not  satisfied,  Pj 
would  decrease  when  n,  is  increased  with  nj  held  constant.  This  should  not  happen  in  an  ideal 
mixture. 

The  calculations  above  should  be  evaluated  with  care.  The  Freundlich  isotherm  does 
not  have  a  proper  Henry's  law  limit  and  is  therefore  not  thermodynamically  consistent.  How 
this  lack  of  a  Henry's  law  limit  affects  the  ultimate  result  cannot  be  determined.  For  this 
reason,  a  similar  set  of  calculations  was  performed  on  a  constant  separation  factor  isotherm. 

Constant  Separation  Factor  Isotherm 

The  Constant  Separation  Factor  (CSF)  isotherm  is  a  mathematically  simple  type  of 
isotherm  of  which  the  Langmuir  is  a  sub-set.  The  CSF  takes  the  form: 


REF'  /  r>  1  \ 

n,  +n,(i?,-l) 


(13) 


The  CSF  has  a  proper  Henry's  law  limit  and  is  thermodynamically  consistent.  As  with  the 
Freundlich  isotherm,  the  CSF  can  represent  favorable  and  unfavorable  isotherm  types  by 
changing  the  value  of  Rj  (Rj  <  1  favorable,  Rj  >  1  unfavorable).  Since  no  analytical  lAST 
solution  was  found  using  the  CSF,  the  lAST  equations  were  solved  numerically.  The  results 
of  the  numerical  solutions  can  be  found  in  Figure  1.  This  plot  shows  how  P2  varies  with  n^ 
while  holding  n2  constant  for  a  particular  set  of  parameters.  P2  is  normalized  to  P2”®®  which 
is  the  pure  component  vapor  pressure  calculated  using  n2  =  0.1  in  Equation  13.  Thus,  when 
is  greater  than  1,  the  theory  predicts  that  the  mixture  partial  pressure  for  component 
2  is  greater  than  that  found  in  pure  component  adsorption  for  the  same  amount  adsorbed  (02  = 
0.1).  Conversely,  when  is  less  than  1,  the  theory  predicts  that  the  mixture  parti^ 

pressure  for  component  2  is  less  than  the  analogous  pure  component  value.  In  all  cases,  the 
amount  adsorbed  for  component  2  (n2)  is  fixed  with  the  other  variables  (nj,  P„  and  P2)  allowed 
to  change.  Since  the  mixture  calculations  performed  assume  ideality,  the  two  adsorbates 
should  exhibit  no  attractive  interactions.  As  a  consequence,  the  only  interaction  the 
components  have  with  each  other  should  be  spatial  interference  or  the  occupation  of  area.  If 
the  surface  is  loaded  with  one  component,  then  the  only  way  to  add  a  second  component  is  to 
displace  an  equal  area  of  the  first  component  to  make  room.  When  the  surface  is  not  fully 
loaded,  the  effect  remains  the  same,  but  is  manifested  in  a  higher  partial  pressure  for  the  first 
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FIGURE  1:  Constant  Separation  Factor  Isotherms  with  Constant  ng 


component  to  "force”  the  original  amount  to  remain  adsorbed.  Therefore,  the  addition  of  a 
second  componrat  (in  this  case  adding  componrat  1)  should  crowd  out  the  original  pure 
component  (component  2)  and  require  a  higher  pressure  to  maintain  nj  at  0.1.  As  the  plot 
shows,  however,  whoiever  one  or  both  of  the  components  demonstrate  unfavorable  adsorption 
(Rj  >  1),  the  mixture  P2  value  is  lower  than  that  found  when  component  2  is  adsorbed  as  a 
pure  component.  Once  again,  it  appears  that  the  two  componrats  are  attracting  each  other  in 
the  adsorbed  phase.  As  with  the  Freundlich  proof  above,  this  should  not  happai  in  an  ideal 
system. 


It  can  be  shown  that,  in  the  limiting  case  where  the  monolayer/saturation  coverages  of 
a  binary  system  are  identical,  lAST  calculations  on  an  unfavorable  CSF  isotherm  result  in  the 
following  equation: 


= 


t  ,  n  „  €,JRT 


'1'  r 


(14) 


As  with  the  v^  similar  form  for  favorable  isotherms  (Equation  2.2),  the  subscripts  can  be 
switched  to  create  the  component  2  equation.  Combining  the  Equation  14  versions  for  both 
components,  performing  a  little  algebra,  and  taking  a  derivative  as  in  Equations  9  and  10 
produces: 


L  = 


^”1  ^  b^e^^{Q+n^+n^ 


(15) 


This  dwivative  must  always  be  less  than  zero.  As  with  the  Freundlich  analysis  above,  values 
for  this  derivative  that  are  less  than  zero  indicate  that  the  theory  is  predicting  an  attractive 
action  between  the  two  components. 

CONCLUSIONS 

The  analysis  of  the  PFCB/water  adsorption  data  indicates  that  it  may  be  inappropriate 
to  apply  the  AST  to  the  specific  adsorption  system.  The  mathematical  analyses  show  that  the 
lAST  should  be  applied  with  care  on  systems  which  demonstrate  unfavorable  single  component 
isotherms.  Indeed,  these  results  sq)pear  to  apply  if  any  component  displays  unfavorable 
adsorption  at  any  point  in  its  single  component  isothom  that  is  integrated  in  the  lAST 
calculations.  This  includes  all  pressures  below  the  mixture  partial  pressure,  and  in  some  cases, 
pressures  above  the  mixture  partial  pressure. 

The  AST  models  mixture  adsorption  by  characterizing  the  adsorbate/adsorbent 
interactions  through  single  component  isotherm  data,  and  providing  adsorbate/adsorbate  (as 
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needed)  interactions  through  solution  thermodynamics.  The  AST  provides  no  mechanism  to 
s^Mirate  adsorbate/adsorb^t  interactions  from  "sdf  interactions"  of  the  adsorbates.  In  the 
case  of  water  adsorption  on  activated  carbon,  or  any  unfavorably  adsorbed  compon^it,  the 
shape  of  the  isotherm  is  a  result  of  a  greater  attraction  between  the  adsorbed  molecules  than 
between  the  adsorbate  and  the  surface.  As  more  is  adsorbed,  the  adsorbate/adsorbate 
attraction  becomes  more  dominant.  Since  the  AST  does  not  separate  this  self  interaction  from 
adsorbate/adsorbent  interactions,  the  self  interactions  will  be  translated  (along  with  all  the 
other  characteristics  which  dictate  the  shape  of  the  single  component  isotherm)  into  the  mature 
equations.  As  a  result,  additional  adsorbates  vdll  reflect  this  adsorbate/adsorb^t  interaction 
and  therefore  demonstrate  an  increased  attraction  to  the  surface  as  more  of  the  unfavorably 
adsorbed  componrat  is  added.  In  cases  where  adsorbate  self  interactions  are  identical  for  each 
component,  and  identical  to  the  unlike  adsorbate/adsorbate  interactions,  this  effect  may  be 
appropriate.  However,  when  componrats  do  not  bdiave  in  this  similar  manner,  tiie  AST 
should  be  applied  with  care. 
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